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ABSTRACT 

The VLB A was used at 1.6 GHz to observe a target field 50' in diameter including the core 
of M31. Novel VLBI correlation techniques were used to observe 200 sources simultaneously, of 
which 16 were detected. We classify all 16 as background AGN based on their X-ray properties 
and arcsecond- and mas-scale morphology. The detected sources were then analyzed for evidence 
of scatter-broadening due to the ionized ISM of M31. The detection of a compact background 
source only 0.25 kpc projected distance from M31 places a constraint on the extent of any extreme 
scattering region associated with center of M31. However, the two sources closest to the core show 
evidence of scatter broadening consistent with that which would be observed if a compact source 
were observed through the inner disk of our Galaxy at the inclination of M31. We interpret 
this as a detection of the ionized ISM of M31 along two lines of sight. With the increases in 
bandwidth and sensitivity envisaged for future long-baseline interferometers, this should prove 
to be a remarkably powerful technique for understanding the ionized ISM in external galaxies. 

Subject headings: Instrumentation: interferometers - Radio continuum: general - Techniques: interfer- 
ometric - Galaxies: individual: M31 - Galaxies: ISM - Galaxies: Active - 



1. Introduction 

Wide-field Very Long Baseline Interferome- 
try (VLBI) is a technique which allows simul- 
taneous observations of multiple sources in a 
single pointing. One of the main applications 
of wid e-field VLBI since its earliest develop- 
ment ( Pedlar et al. 19991) has been the study 
of nearby galaxies , particularly starb urst galax- 



ies such as M82 (IFenech et al. 2010), NGC253 



( Lenc fc Tingavl2006l ). andNGC4945 jLenc fc Tingav 
20091 ). Such studies have provided insights into a 
wealth of astrophysical processes, providing esti- 
mates of the supernova rate in these galaxies, and 
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the electron density of the ionized medium that is 
giving rise to the starbust. 

M31 has not, until now, been the subject of such 
a study; we ascribe two reasons for this. Firstly, 
the large angular extent of the galaxy means that 
it is only with recent developments in VLBI corre 



lation and data reduction techniques (jPeller et al 



20111: Im organ et al.l 120111) that it has been possi- 
ble to survey any more than the tiniest fraction 
of M31 using VLBI. Secondly it is not a starburst 
galaxy and so is unlikely to have the large popu- 
lation of young Supernova Remnants (SNRs) and 
exotic objects present in these more extreme envi- 
ronments. 

However, with the increasing sensitivity of 
VLBI due to bandwidth upgrades such as that cur- 
rently underway at the VLBA0, compact sources 
intrinsic to M31 may be observable. For exam- 
ple at the distance of M31, the Pulsar Wind 



'Very Long Baseline Array 
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Nebula (PWN, iGaensler fc S J l2006l) associ- 
ated with the Crab pulsar would have a total 
flux of 6mJy and an extent of approximately 1" 
with filamentary struc ture on a scale of 15 mas 
(jBietenholz et al.lll997() , well matched to the res- 
olution and sensitivity of 1-2 GHz VLBA obser- 
vations. There are also objects similar to PWNe 
but which are powered by microquasars, for ex- 
ample the W50 ne bula associated with SS443 
( Dubner et alllQQSl) . An object such as this would 
also be detectable in our observations. In addition, 
an ultra-luminous stellar-mass mi croquasar within 



our ta rget field was reported by iMiddleton et al 
(I2ni3[ ). Wide-field VLBI studies of M31 therefore 



have potential to detect analogs of rare Galactic 
objects in other galaxies. 

VLBI measurements of the proper motion of 
M31 are currently lacking, in co ntrast to many 



other nearby gal axies such as M33 ([Brunthaler et al 



20051 ) and ICIO (|Brunthaler et al.ll2007l) . Deriving 
such proper motions via VLBI techniques requires 
both compact sources intrinsic to the galaxy (usu- 
ally, though not necessarily, masers), and back- 
ground compact sources which can be used as 
calibrators. Discovery of the latter is yet to be re- 
ported. Such proper motion measurements would 
be a useful complement to the recent proper mo- 
tion for several fields observe d at multiple epoc hs 
with the HST as reported bv lSohn et"al] (|2012l ). 



1.1. M31 

M31 is the nearest non-dwarf n eighbor to the 
Milky Way, at a distance of 785 kpc (jMcConnachie et al, 
20051 ). Its D25 isophote is an ellipse with major 
and minor axes of 155' and 50' respectively, and 
the inclination of the disk to our lin e of sight 
is approximately 77° (|Tu11v fc Fisheil [XOSSI. In 
all continuum surveys and in far-infrared images, 
diffuse emission from a central region some 4' (ap- 
proximately 1 kpc) in extent is seen. The other 
notable feature both in low-resolution radio maps 
and in the far infrared is a ring some 10 kpc dis- 
tant from the center (see Figure [T]) . A rather less 
obvious but nonetheless important feature of M31 
was revealed in Ha , N li and O i l l ima ging of the 
core carried out bv lJacobv et al.l ( 19851 ). They re- 
vealed a striking spiral structure (with an entirely 
different orientation to the main galactic plane) in 
the ionized gas occupying this central region. 

M31* is the M31 analog of Sgr A*: an X-ray 
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Fi ^. 1. — Gra y scale: A subset of data presented 
Beck et all (|l998l) . extracted from the VLA 



archive. The core and 10 kpc ring are clearly vis- 
ible. Contours: predicted RMS sensitivity of our 
observation (|J.Jy). The peak sensitivity is 30|J^Jy 
RMS 



and radio source corresponding to a supermassive 



black hole of mass 1.4;°] X IO^Mq (iBender et al 



20051) (c.f. Sgr A* - 4.0 x IO^Mq). It is three 
times less luminous than Sgr A* in the radio and 



three times more luminous in X-rays (jCarcia et al 



201fl[ ). significantly below the fundamental plane 
The estimated star formation rate in M31 is sig- 
nificantly less t han in nearby starburst galaxies. 
Williama ( 2003 ) found the optical mean star for- 
mation rate for the disk of M31 to be ~1 M© yr~^, 
approximately half that found in M82 from radio 
observations; while iKang et al. (|2009l ) found an 
average of 0.6 to 0.7 Mq yr-i over the last 400 Myr 
from UV observations. One supernova has been 
obser ved in M31, very close to the center ( Hartwi j 
18851 ): and a number of SNRs in the co re have been 



identified via their X-ray properties (jKong et al 



20031 ) ■ and throughout the North-Eastern half of 
the galaxy via their forbid den line characteristics 
( Braun fc WalterbodHigl . 

Unsurprisingly, there exists a great deal of 
archival data on M31. Radio continuum surveys 
are summarized in Table [Tl The h i ghest resolution 
survey to date is that of 



Braun (|l990l . hereafter 
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B90) which imaged just over one half of M31 using 
10 pointings of the VLi^ at 1 .4 GHz with a resolu- 
tion of 5" using the B-, C-, and D-configurations. 
This survey detailed 534 individual continuum 
sources, a large fraction of which were unresolved 
and of which 103 have possible associations with 
Ha nebulae and seven with previously identified 
SNRs. 

The resulting map also contains 121 objects as- 
sociated with previously known radio sources de- 
tected in surveys of t he entire galaxy at 



610 MHz (iBvstedt et al.M l984D and at 1412 MHz 



(jWalterbos et al.Ml985[ ). the 3 6W and 37W cat 



alogs respectively. In addition, iBeck et al.l (|1998l ) 
surveyed the entire galaxy at 20 cm using both the 
Effelsberg dish and the VLA in the compact D 
configuration, concentrating on polarization prop- 
erties of the M31 field. Su bsequently the entire 
M31 field was surveyed by iGelfand et all (|2004l ) 
(hereafter GLG) with a single pointing of the VLA 
in A-configuration at 325 MHz, uncovering a total 
of 405 sources across 7.6 deg^. 



In the X-ray, IStiele et al.l (|2011l ) used archival 
data, supplemented with new observations to pro- 
duce a deep XMM-Newton Survey of the entire 
disk, producing a catalog of 1987 sources. 

In H I emission stu dies, the stat e -of-th e-art is 
the extensive study of iBraun et al. (|2009[ ). Sur- 
ve ys of H ii regions have bee n un dertaken both 



by IWalterbos fc BraunI (|1992| ) and lAzimlu et al 
(l201lh . 



We have undertaken a pilot study to explore a 
field of M31 with wide-field VLBI, in detail, for 
the first time. In this paper we present these ob- 
servations, supplemented by some data from the 
VLBA archive and VLA observations. In Section[2] 
we describe the observations and our source selec- 
tion and data reduction procedure; in Section[3]we 
describe the results of our observations and cross- 
match our sources with previous surveys and data 
from the VLBA archive; in Section 2] we discuss 
our results. 



■^Very Large Array (now upgraded and renamed the Karl G. 

Jansky Very Large Array) 
^The Westerbork Synthesis Radio Telescope 



2. Observations 

An 8-hour observation was carried out with 
aU 10 antennas of the VLBA on 2010-July-04. 
The target field was centered on 00'^42'"43.3% 
41°28'25" (J2000.0). It was chosen using the fol- 
lowing criteria: 1) overlap with the B90 survey; 
2) overlap with both the core and a segment of 
the starforming ring; 3) inclusion of bright, un- 
resolved sources which provide a high probability 
of easy detection (Figure [1]). The source ICRF 
J003824. 8-1-413706 was used as a phase reference 
calibrator, with a duty cycle of 5 minutes on tar- 
get, 1 minute on calibrator. The source 3C84 was 
used as a fringe finder and bandpass calibrator and 
was observed for 5 minute scans 3 times during the 
observation. 

64 MHz of contiguous bandwidth between 
1608.49 and 1672.49MHz was recorded in 4x16 
MHz sub-bands. Both polarizations were recorded 
with 2-bit precision for a total data rate of 
512 Mbps per antenna. 

2.1. Target Selection and Correlation 

Correlation was c arried out u s ing t he DiFX 



software correlator (jPeller et al.l 120071). The 



multi-phase center technique (jPeller et al.l 120111 ) 



was vital for correlation of the target field as this 
allowed multiple phase centers to be correlated 
simultaneously in a single pass. 

The remarkable efficiency of the DiFX correla- 
tor in multi-phase center mode means that there 
is little processing time to be saved in being eco- 
nomical with the number of phase centers. Con- 
sequently, we did not consider whether a source 
showed a compact morphology, or if it was bright 
enough to be detected (while the sensitivity of our 
observations equal that of the B90 survey at our 
pointing center, primary beam effects mean that 
this is certainly not the case over the entire field). 
This was justified on the basis that any source 
might be variable, and any resolved source could 
have a compact component and/or be masking a 
separate compact source. 

All B90 sources closer than 20' to the field cen- 
ter were targeted. Use was also made of the 37W 
catalog, which covers a wider area and also gives 
separate coordinates for each source component: 
all 37W source components within 30' were tar- 
geted. In addition, sensitive VLA images of the 
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Table 1 

Radio Continuum Surveys of M31. 



Survey 


Telescope 


Frequency 


Resolution 


RMS Sensitivity 


Reference 






(GHz) 


(") 


(mJy) 








B90 


VLA-B,C,D 


1.5 


5 


0.03 


uniform 


Braun (1990) 




37W 


WSRT 


1.4 


23x35 


0.2 


peak 


Walterbos et al. (1985) 


36W 


WSRT 


0.61 


54x82 


0.8 


peak 


Bvstcdt et al. (1984) 


GLG 


VLA-A 


0.325 


6 


0.6 


peak 


Gclfand et al. 


20041 


NVSS 


VLA-D 


1.4 


40 


0.45 


uniform 


Condon et al. 


T998) 


WENSS 


WSRT 


0.325 


54x82 


3.6 


uniform 


R,engelink et a 


. (1997) 


VLSS 


VLA-A 


0.070 


80 


100 


uniform 


Cohen et al. (2007) 



Note. — All but B90 cover the entire disk. 



core reg ion of M31 were extr acted from the NVAS 
archive (jCrosslev et al.ll2007l ) , however no firm de- 
tections could be discerned other than those al- 
ready included in the B90 catalog. Any source 
which was within 10" of another was removed 
from the list of correlation targets. The output of 
the correlator consisted of 195 separate visibility 
datasets, each phase-centered on one of our target 
sources. 

2.2. VLBI Data Calibration 

For the most part, calibration of a wide- field 
VLBI dataset proceeds in an identical fashion to a 
standard VLBI observation. First, extensive flag- 
ging was carried out via visual inspections of the 
visibilities. Next, corrections to the amplitude 
were made using th e AIPS task 'ACCO R'. Our 



own ParselTongue ( Kettenis et al. 20061 ) scripts 



were used to apply phase corrections based on 
pulse cal tone information extracted by DiFX. 
DiFX extracts all pulse cal tones from across the 
band, rather than just the edge tones, increas- 
ing the distance between 27r ambiguities of phase 
jumps to l[xs. Visual inspection of the visibili- 
ties confirmed that this ensured that all sub-bands 
were coherent with each other throughout the ob- 
servation (in contrast to the standard VLBA pro- 
cedures, which were unable to correctly adjust for 
phase jumps of > 0.07ns). 

Next, system temperature measurements and 
gain curves were used to provide amplitude cal- 
ibration using the AIPS tasks 'APCAL'. Online 
system temperature measurements for the Han- 
cock antenna were found to be invalid and so a 
nominal value of 40 K was used throughout. 



A combined fringe fit solution for all sub-bands 
was then derived using the AIPS task 'FRING', 
assuming a point source calibrator and a fringe 
rate of zero. The resulting calibration solution was 
then applied for a first imaging attempt. 

2.3. Wide-field Data Calibration and Imag- 
ing 

Before applying any calibration, phase centers 
for the 'duplicate' sources removed from the corre- 
lation target list (see Section [2T|) were regenerated 
using the AIPS task 'UVFIX' to shift the phase 
center of the nearest correlated phase center to 
this point. This resulted in 226 visibility datasets 
consisting of 147 B90 sources and 86 components 
of 79 sources from the 37W catalog. 

The calibration solution as derived in the pre- 
vious section was then applied to all these sources. 
In addition primary beam c orrections as described 
bv iMiddelberg et ah ( 2013 ) were applied. In this 
scheme, the primary beam is assumed not to vary 
across the area which is to be imaged, allowing 
the primary beam corrections to be applied di- 
rectly to the visibilities, taking into account the 
slight (but significant) dependence of the correc- 
tion on antenna, polarization, and hour angle. 
The frequency scaling of the primary beam was 
calculated for each s ub-band. A primary b eam 
model as described bv lUson fc Cotton ( 20081) was 
used with parameters measured for each individ- 
ual VLB A antenna (R.C. Walker Priv. Comm.; 



see also IMiddelberg et al.| Table 3). Suc h cor- 



rections were verified bv Middelberg et al. to be 



accurate at a level of 4%. 
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Preliminary images of all target sources were 
then produced. The RMS noise level in the im- 
age was determined using the AIPS task 'IMEAN'. 
Comparison of this value to the image peak value 
was sufficient to easily identify the brighter of the 
detected target sources. 

2.4. Refinement of the Calibration 

For both datasets the coherence time was far 
greater than the phase reference duty cycle for the 
vast majority of the observation, given the close 
proximity of the calibrator to the target fields and 
the favorable time of the observation in the solar 
cycle. 

Further careful flagging was carried out, in- 
cluding flagging all data whose u coordinate was 
< 50 kA. Such visibilities are susceptible to RFI 
due to the low fringe rate and this flagging was 
successful in removing horizontal stripes that had 
been seen in some images. 

The detection of bright sources within the 
target field allows the possibility of using these 
sources as in-beam calibrators. With sufficient sig- 
nal to noise (S/N), these sources can be expected 
to give improved calibration solutions due both 
to their proximity to the other targets and the 
fact that they are observed simultaneously rather 
than in interleaved pointings. One option was to 
use the multi-s ource phase referencing technique. 



as described in iMiddelberg et al.l (j2013[ ). however 



the S/N of all our sources combined in quadrature 
would only be approximately 10% greater than 
for the brightest source alone, which was deemed 
unlikely to lead to a significant improvement. 

The use of 37W144 as an in-beam calibrator 
(phase only, with one combined solution for all 
sub-bands and polarizations) resulted in an aver- 
age increase of S/N of most of the target sources 
of 10%. 

This refined calibration was then applied to all 
target visibility datasets and the data were then 
averaged in frequency to 8 spectral channels per 
sub-band. These averaged data were used for all 
subsequent analysis. Source detection images were 
then produced for all fields consisting of 'dirty', 
normally- weighted maps. 



2.5. Source Detection 

Since the difference in synthesized beam size 
between the B90 observations and our VLBA ob- 
servations is approximately 1000, a minimum of 
1 million resolution units must be searched for each 
source. Imaging requirements meant that in real- 
ity 8192x8192 pixel images with each pixel mea- 
suring 1 mas on a side were used to image the 
5" resolution unit of the B90 catalog. Assum- 
ing a Gaussian distribution of pixel brightness, we 
would expect one 5.6 cr noise peak in every VLBI 
image. We are therefore motivated to determine 
the probability that any apparent source is real, 
rather than just a noise spike. 

We calculate a likelihood ratio (LR) for each 
pixel in each of our images formed from the fol- 
lowing hypotheses: 1) the observed flux can be 
associated with a cataloged source (i?i); 2) the 
observed flux is a noise spike [Hq). This calcula- 
tion is made under the following assumptions: 1) 
the noise in each image is Gaussian-distributed; 
2) the RMS of the noise does not vary across any 
individual VLBI image. 

Hf) is simply the Cumulative Distribution Func- 
tion of the positive half of the Gaussian distribu- 
tion: 



p{m,arms]Ho) = CDFA^(m, CTrms) 



1 - erf 



2 

rms . 



(1) 

(2) 



where m is the brightness of the pixel in the VLBI 
image and arms is the RMS noise in the VLBI 
image. 

Hi is based solely on the position of the VLBI 
pixel and the cataloged source position/morphology. 



p{r,d;Hi) = l{r) 

= (^(r,CTa,Crf,,( 



(3) 
(4) 



where r is the relative vector offset between cata- 
loged position and VLBI pixel position. The cat- 
aloged source structure 6 is described by aa, <yb 
and where these are the cataloged semi-major 
axis, semi-minor axis and position angle of the 
source respectively. l{r) is then calculated using a 
normalized elliptical Gaussian, iy9, parameterized 
by these variables. The instrumental resolution is 
also accounted for here as detailed below. 
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The Likelihood ratio (LR) is then simply 
P{Hi) 



LR: 



(5) 



A further refinement we considered was taking 
into account how well-matched the cataloged flux 
and the VLBI flux were. For example, the cata- 
log/VLBI flux scaling could have been modeled as 
a top-hat function convolved with a softening that 
incorporated image noise and variability softening. 
This was not done for several reasons. Firstly, 
any source could be variable — particularly on the 
20-year timescale that separates our observations 
from the catalog observations. Secondly, we ex- 
pect that a significant fraction of the flux could 
be resolved out by the VLBA. This leaves only 
cases where the VLBI detection is significantly 
brighter than its cataloged counterpart. Since the 
object of this exercise is to discriminate between 
noise spikes and weaker sources, such a refinement 
would serve no purpose. 

2.5.1. Generating LR Maps 

While the quoted resolution for the B90 obser- 
vations was 5" the astrometric precision could be 
expected to be greater for brighter sources. For all 
B90 sources, cTq and ct;, were set either to their cor- 
responding values, or to the quoted error on these 
values, whichever was larger. 

Likewise, IWalterbos et al. I (|l985l) gave a major 
axis, minor axis, and a positional angle </> for each 
resolved source and an error in RA and declination 
for all sources. The maximum values were used for 
a a and Cf,. 

FITS images corresponding to the VLBI images 
containing the LR values for each pixel were then 
generated. Before any use was made of the VLBI 
maps, the outer 512 pixels from the edge were ex- 
cluded from the analysis to mitigate the effect of 
aliasing effects. 

2.5.2. Determining the Detection Threshold 

Table [2] gives the VLBI map statistics along 
with maximum value from the LR maps for all 
detected sources and the field containing the next 
highest LR pixel. 

We note that the simple statistic of VLBI S /N 
is as good a predictor of whether a source is real 
as the likelihood ratio. This reflects the extremely 



poor constraint on position provided by the cata- 
logs, compared with the VLBI position. 

These data show also that our assumption of 
pure Gaussian noise is a good one. Our brightest 
noise spike in ^ 170 x 7168 x 7168 pixels is 6.15 a. 
A pixel of this value or brighter would be expected 
to occur in Gaussian noise once in every 50 maps. 
Compare this with our weakest detection which 
would only occur in Gaussian noise once in every 
64 000 maps. It should be stressed that such pure 
Gaussian noise in the non-detection maps was only 
achievable with very careful flagging of the visibil- 
ity data and with masking of the VLBI images to 
mitigate aliasing effects at the image edge. 

This very clear cutoff between our weakest de- 
tection and our brightest noise peak was conflrmed 
upon visual inspection of the VLBI images and 
LR maps. We note, however, that real sources, 
weaker than 6.15 cr almost certainly exist in our 
data, particularly in those flelds near the edge of 
the primary beam, however we flnd no statistical 
evidence for them in this analysis. 

For the upper limit on the brightness of any 
non-detection, we adopt 6.6 cr. This corresponds 
approximately to P< 0.05 of a source being a noise 
spike over ^ 200 images. 

2.6. Wide-field Mapping 

In addition, each targeted source was imaged 
with visibilities on baselines longer than 10 MA dis- 
carded (approximately corresponding to the inner 
8/45 baselines). This allowed the pixel size to be 
increased to 5 mas with a commensurate increase 
in field size to a 40" diameter. The purpose of 
this exercise was twofold: first, to detect any emis- 
sion from extended parts of the sources; second, to 
search for less compact sources which would not be 
detected on the longer baselines. 

The same likelihood analysis was then applied 
to these images. No new sources were detected. 
Though the S/N of some of the more compact 
sources dropped, the sources with the highest LR 
values were the 16 sources already detected (see 
Table O. 

2.7. Imaging and UV Modeling of De- 
tected Sources 

After source identification, new calibrated visi- 
bility datasets were generated, centered on the po- 



6 



Table 2 
Detection Statistics 



All baselines 10 MA cutoff 



"ipcak tT,™ S/N -logio(p(^fo)) logio(LR) S/N -logio (p(ffo)) log^o (LR) 

(mjy) (/jJy) (mpcak/o-rmn) 



11.95 


65 


180 


>308 


4.64 


59 


79 


>308 


4.70 


79 


60 


>308 


3.14 


64 


49 


>308 


0.73 


34 


21 


102 


1.64 


80 


20 


92.6 


2.35 


119 


20 


85.7 


1.96 


102 


19 


81.6 


17.46 


1074 


16 


58.9 


0.75 


58 


13 


37.6 


0.68 


53 


13 


37.3 


2.68 


272 


9.8 


22.4 


0.50 


53 


9.5 


20.8 


0.38 


42 


9.0 


19.1 


0.53 


62 


8.5 


17.2 


0.40 


55 


7.2 


12.5 


0.53 


87 


6.2 


9.41 



>308 


240 


>308 


>308 


>308 


64 


>308 


>308 


>308 


77 


>308 


>308 


>308 


41 


>308 


>308 


102 


22 


111 


110 


92.0 


15 


50.0 


49.3 


85.2 


13 


39.1 


38.5 


80.7 


14 


45.1 


44.1 


56.1 


10 


26.0 


23.2 


37.0 


8.5 


16.9 


16.4 


36.8 


14 


46.4 


45.9 


22.0 


8.0 


15.3 


14.9 


20.3 


6.5 


10.3 


9.9 


18.7 


6.9 


11.5 


11.1 


16.7 


9.6 


21.6 


21.1 


11.4 


8.8 


18.2 


17.7 


8.94 


6.3 


9.81 


9.37 



Note. — The 16 firm deteetions are shown, along with the non-detection with the highest LR, (below the 
horizontal line). These are different fields for the two different resolutions, mpcak, p(^o): s,nd LR are shown 
for the brightest pixel in the image. 



sition of the brightest pixel in the LR maps. The 
package DIFMAP was then used to fit a circular 
Gaussian model to the visibilities. The visibilities 
were finally re-shifted to the centroid of the Gaus- 
sian fit and final images were produced (Figure [2]). 
Flux measurements and their uncertainties 



were then derived using BLDBCAT (IHales et al 
I2OI2I) with the same parameters used bv lMiddelberg 
(j2013l ) (i.e. a pixellation error of only 1% since 
the pixels are much smaller than the beamsize, 
and assuming a nominal surface brightness er- 
ror of 10%). No errors associated with the pri- 
mary beam correction are explicitly accounted 
for. 37W060 is ju s t bey ond the radius probed by 
iMiddelb erg et al.l (l2013l ). and given the primary 
beam correction factor of > 23, it is possible that 
the flux error is greater than the nominal 10%. 

To characterise the source morphology, it was 
necessary to produce a robust estimate of the 
source size. These visibilities were then averaged 
in frequency and polarization and averaged in time 
over 2 hours to produce approximately 200 visibil- 
ities per source. These datasets were used to fit a 
Gaussian to the real part of the visibilities against 
uv distance (corresponding to a circular Gaussian 



fit in the image plane). The visibility weights and 
the scatter in the imaginary components were used 
to estimate the error bounds of each individual vis- 
ibility. 

The error on the fit (both fiux and size) was 
then estimated using a perturbation analysis: a 
random amount of noise was added to the real 
et al.l component of each visibility (drawn from a Gaus- 
sian distribution with the standard deviation given 
by the error bound of that visibility) and the fit 
was recalculated. This process was repeated 1000 
times and the standard deviation on the distribu- 
tions of the fiux and the source width in the image 
plane provided the errors on these quantities. 

The final source positions and their errors are 
shown along with the BLOBCAT fiux, fitted Gaus- 
sian size and zero-spacing flux in Table [3l The fit 
to the visibilities is shown in Figured 

2.8. Source Astrometry 

The astrometric error on the source positions 
contains two components. The first component 
is the formal error given by the limited S/N of 
the Gaussian fit. The formal error for the posi- 
tion component of a Gaussian fit is given approx- 
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imately by 



fir 



(6) 



"pos — >^pos 2^ 

(see lReid et"allll988f) where Spos was the FWHM 



of the Gaussian fit, or the beamsize, whichever 
was largeiQ, S was the zero-spacing flux from the 
Gaussian fit and (Trms was the RMS noise in the 
source-detection map. 

However, this formal error considers only the 
accuracy to which a Gaussian model can be fitted 
to the data and ignores any systematic shifts in 
the data itself. For phase referenced VLBI obser- 
vations at 18 cm, there are two primary effects: 

1. The positional error of the primary calibra- 
tor source; and 

2. Large-scale structure in the ionosphere at 
the time of the observation, which leads to 
phase referencing errors, which displace the 
target sources. 

The primary phase reference source, ICRF 
J003824. 8-1-413706, has a catalog positional accu- 
racy of 0.1 mas. However, the ICRF positions are 
defined by observations at 2.3 and 8.4 GHz, and 
are known to be affected in some cases by source 
struc ture, both tim e invariant and time variable 
(e.g.,|Porcai|200i). As shown in Figure El ICRF 
J003824.8-H413706 clearly has a core-jet morphol- 
ogy and is likely to be affected by this issue, espe- 
cially since our observations are made at the lower 
frequency of 1.6 GHz. Therefore, we consider it 
likely that the reference position used for a point- 
source model of ICRF J003824.8-t-413706 at 1.6 
GHz is in error by 1-2 mas. 

The magnitude of phase referencing errors is 
difficult to estimate, since they depend on the (dif- 
ficult to model) ionospheric conditions at the time 
of the observation. However, since the target- 
calibrator separation is quite small, only ~1°, this 
contribution is likely to be relatively small. Based 
on astrometric observations of pulsars at 1.6 GHz, 
a value of 0.25 - 0.5 mas mi ght be expected for a 



given observing epoch (e.g. , iBrisken et al.l 120021 ) 



Thus, the total systematic error contribution is 
likely to be of order 1-2 mas, dominated by the 



*The beamsize for our observations was 9.77x5.55 mas with 
a PA 0.37° therefore 5.55 mas was used as the minimum 
5pos for the Rj 
hnation error. 



uncertainty in the position and structure of the 
reference source ICRF J003824.8-I-413706. 

For all sources, the formal position error esti- 
mate is smaller than this systematic error esti- 
mate, by a large margin in the case of brighter 
sources. The use of the source 37W144 as an in- 
beam calibrator does not change these systematic 
error contributions, since the position derived for 
37W144 was subject to these errors in the ini- 
tial imaging. In essence, these errors determine 
the offset of 37W144 from its true position, and 
the subsequent use of 37W144 as an in-beam cal- 
ibrator then reinforces this offset. However, when 
considered relative to 37W144, the residual er- 
rors in positions for the other sources are much 
smaller. Since the angular displacements are <20', 
and there is no time interpolation of the solutions, 
the residual systematic errors are probably <0.1 
mas, again drawing on ex perience from pulsar as - 
trometry at 1.6 GHz (e.g., Chatteriee et al.|[2009() . 
For many sources, this is smaller than the formal 
fit uncertainty. 

Accordingly, in Table |3] we list only the formal 
fit errors for positional quantities, and note that an 
additional uncertainty of 2 mas should be added 
in quadrature to each axis to obtain an abso- 
lute position error. If future observations improve 
the accuracy of the 1.6 GHz position of ICRF 
J003824.8-K413706 (and by extension 37W144), 
this error component could be reduced. 

2.9. Reduction of archival VLBA data 

The VLBA archive contains data from an 8- 
hour phase- referenced observation of 37W142. 
The observation took place on 2003-09-19 and the 
same phase reference source was used as for our 
observations. In addition 37W144 was observed 
intermittently as a check source. Four 4 MHz sub- 
bands were recorded, two at 18 cm and two at 
21 cm. Data were available for 9 antennas (Brew- 
ster being the exception). 

Data reduction proceeded in much the same 
fashion as for our main observation with a cou- 
ple of exceptions: before calibration, earth orien- 
tation parameters were updated with the correct 
values. Other than the removal of all Owens Valley 
baselines, no flagging was deemed necessary and 
standard VLBA pulse cal corrections and system 
temperatures were used for phase reference cali- 
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Fig. 2.— Normally-weighted images of the 16 detections. In order, (a) 37W060, (b) 37W080, (c) 37W092, 
(d) 37W093, (e) 37W115, (f) 37W118, (g) 37W123, (h) 37W125, (i) 37W129, (j) 37W142, (k) 37W144, (1) 
37W150, (m) 37W156, (n) 37W158, (o) B23, (p) B157. The first contours are at ±3-sigma, increasing in 
increments of ^f2. 
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Fig. 3.— Plots of the model fit of a Gaussian to the visibilities. In order, (a) 37W060, (b) 37W080, (c) 
37W092, (d) 37W093, (e) 37W115, (f) 37W118, (g) 37W123, (h) 37W125, (i) 37W129, (j) 37W142, (k) 
37W144, (1) 37W150, (m) 37W156, (n) 37W158, (o) B23, (p) B157. The red dots show the real parts of the 
visibilities and their associated errors, the blue dots show the imaginary parts of the visibilities (with the 
same errors). The solid black line shows the least squares fit. The two dotted lines show this same fit with 
the standard error on both the zero-spacing flux and the Gaussian width first added, then subtracted. 
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bration. After fringe fitting, the 18 cm and 21cm 
data were separated. Several rounds of phase and 
amphtude and phase self calibration were carried 
out on the phase reference source and this was 
applied to the other sources. The correlation po- 
sition of the calibrator differed by approximately 
0.3 mas from that of the 2010 observations. This 
position was assumed correct for this epoch. 

37W144 was clearly detected and the steps de- 
scribed in Section 12.71 were carried out. 37W142 
was marginally detected (which is consistent with 
our results) although the model fit failed to con- 
verge. 

The results are shown in Table U) 
3. Results 

3.1. Cross-matching virith other surveys 

Counterparts to the detecte d sources were 
searc hed for in the B90, 37W (IWalterbos et al. 
19851). GLG (i Gelfand et all 1200^ and Stiele 



(jStiele et al. 20nJ catalogs. 



The VLBI positions are compared with the po- 
sitions of the counterparts in four surveys in Fig- 
ure m The B90 and the 37W catalogs appear to 
have systematic errors in the position in excess of 
the listed errors in RA and declination. This is 
unsurprising given the considerable difficulties in 
achieving high ast rometric accuracy in wide-field 
imaging (see, e.g. ( Condon et al. 19981) . sect. 4.3). 
However, they remain small relative to the mini- 
mum size we used for our source-finding analysis. 

Using our source detection method with an 
input catalog with high astrometric accuracy, it 
would be possible to detect sources significantly 
below our detection threshold of 6.6a. 

3.2. Source Classification 

Deciding whether the detected sources are in- 
trinsic to M31 or background sources is vital if we 
are to draw any conclusions from our observations. 
The only persistent background sources with 
brightness temperatures in excess of 10^ K (and 
thus easily detectable by VLBI) are AGN. Detect- 
ing 16 sources within on e pointing of the VLBA i s 
broadly consistent with iMiddelberg et all (|201l[ ). 
so it would not be surprising if all of the detected 
sources were background AGN. Indeed, all avail- 
able data are consistent with this being the case. 



In this section we consider the evidence for these 
sources being either intrinsic to M31, or back- 
ground AGN. 

All source s with X-ray count erparts have been 
classified by IStiele et all (|201ll) as AGN on the 
basis of their X-ray properties, on which basis we 
rule out the possibility that any of our detected 
sources are X-ray binaries or associated nebulae. 

Young Supernova Remnants SNRs undergo free 
expansion for some time before beginning to de- 
celerate. Although there is much variation in 
measured expansion velocities, a typical expan- 
sion speed is 10 000 km s^^. This corresponds to 
~0.5mas year~^ at the distance of M31. Thus 
even the very youngest SNRs that are nonetheless 
old enough to be cataloged in B90 would likely be 
resolved into a ring or partial ring by the VLBA. 
In addition, such recent supernovae would very 
likely have been observed. Furthermore, SNRs 
older than ^100 years would begin to be resolved 
out completely, though spots of compact structure 
within the ring could remain detectable on VLBI 
scales. 

Pulsar Wind Nebulae (PWN) are another class 
of Galactic objects to which our observation would 
be sensitive. As already discussed, the Crab neb- 
ula at the distance of M31 would have a total flux 
of approximately 6mJy and a total extent of 1" 
with filamentary structure on a scale of 15 mas, as 
would a nebula associated with an X-ray binary. 

However, none of our sources have a ring or 
partial ring morphology, and with one excep- 
tion (37W115, which has AGN morphology), 
none show extended structure on scales greater 
than ^10 mas. Moreover, with the exception of 
37W115, the VLBI flux matches those from low- 
resolution surveys to within a factor of two, mean- 
ing that little or no extended structure is resolved 
out on scales to which the VLBA is insensitive 
(> 1"). While PWNe vary greatly in morphology 
and radio luminosity, this is smaller than the Crab 
nebula by two orders of magnitude. 

While we cannot entirely rule out the possibil- 
ity that one or more of our detections is a source 
intrinsic to M31, perhaps of a source population 
which does not exist in our own Galaxy, there is 
no real evidence for this. Furthermore, many of 
the detected sources have features which are either 
suggestive of AGN, or which are incompatible with 
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Table 3 

Detected Source properties and Cross-match with Other Surveys 



Name 


RA (J2000.0) 


Dccl. (J2000.0) 


Iblob 


Imodol 


Model Size 


Braun 


37W 


GLG 


Stiele 










(mJy) 


(mJy) 


(mas) 










37W144 


00:42:51.6222130(12) 


41 


:26:34. 318347(23) 


20.6±1.0 


13.20±0.12 


2.67±0.11 


103 


144 


Oil 


1055 


37W123 


00:42:25.919929(21) 


41 


:29:18. 46241(42) 


1.761±0.097 


1.26±0.15 


9.02±2.1 


69 


123 






37W156 


00:43:18.635967(21) 


41 


:24:22. 82067(42) 


0.397±0.058 


0.38±0.07 


0.51±2.2 


131 


156 






B157 


00:43:38.685540(22) 


41 


:30:26. 56947(43) 


0.649±0.073 


0.53±0.10 


2.75±2.6 


157 




025 


1272 


37W125 


00:42:26.656261(16) 


41 


:18:04. 94893(18) 


1.70±0.11 


1.13±0.17 


8.63±2.1 


70 


125 


002 




37W150 


00:43:04.037813(31) 


41 


:38:46. 65434(35) 


1.44±0.10 


1.46±0.52 


19.31±13 


118 


150 




1123 


37W092 


00:41:45.1683974(28) 


41 


:26:23. 584359(56) 


7.63±0.39 


4.94±0.11 


2.20±0.29 


25 


092 


021 


706 


B23 


00:41:41.810643(19) 


41 


:29:36. 45939(38) 


0.896±0.086 


0.84±0.08 


0.00±0.40 


23 








37W129 


00:42:32.252269(11) 


41 


:42:36. 37584(21) 


2.67±0.17 


1.76±0.15 


2.27±1.3 


71 


129 


046 




37W142 


00:42:50.383045(25) 


41 


:16:10. 13693(28) 


2.77±0.16 


1.44±0.30 


14.35±4.1 


100 


142 


001 




37W158 


00:43:26.3158109(39) 


41 


:19:11. 787621(77) 


5.09±0.27 


3.90±0.13 


4.17±0.36 


140 


158B 


005 


1224 


37W060 


00:40:35.778378(13) 


41 


:35:10. 55738(26) 


21.4±1.7 


19.22±1.5 


1.05±0.92 




060 


069 




37W080 


00:41:18.052819(12) 


41 


:23:17. 30580(25) 


2.80±0.20 


2.24±0.16 


0.00±0.79 


4 


080 






37W093 


00:41:46.146853(11) 


41 


:16:39. 24342(22) 


3.15±0.20 


2.28±0.20 


3.47±1.0 


29 


093 


016 


715? 


37W115 


00:42:18.7006175(37) 


41 


:29:26. 297058(74) 


7.45±0.39 


5.08±0.18 


2.29±0.46 


61 


115 


019 


856 


37W118 


00:42:21.983656(21) 


41 


:08:00. 99797(36) 


5.37±0.43 


3.66±0.64 


5.85±2.8 


66 


118 


008 





Note. — Ibiob and I„ 
(see Section |3.2.1L 



ci arc the model and BLQBCAT flux as described in Section [2T7] The cross-match with Stiele 715 is only tentative 



Table 4 
Archival Data on 37W144 



Souree 


Wavelength 






Iblob 


Imodel 


Model Size 




(cm) 


(mas) 


(mas) 


(mJy) 


(mJy) 


(mas) 


37W144 


18 


0.80±0.15 


-0.19±0.34 


lO.lil.l 


11.00±0.82 


1.5±1.1 


37W144 


21 


0.77±0.21 


0.03±0.55 


11.1±1.3 


15.3±2.1 


4.2±2.3 



Note. — Sb,a and (5Dec give the difference of the fitted position offset from the positions derived 
from the sources in our own observations. The error is for the fit on the archival data only. 
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Fig. 4. — OflFset in position of sources in (a) Braun, (b) 37W, (c) GLG, (d) Stiele compared with VLBI 
positions. West is positive. The centroid of the gray hnes is the weighted moan of the positional offsets. The 
width of the gray hne is the cataloged errors for all sources combined in quadrature. 37W115 is excluded 
due to its arcsecond-scale structure 
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the galactic sources that we might expect to detect 
in M31. These are summarized for each individ- 
ual source below within four general categories: 
background AGN on the basis of large-scale mor- 
phology; background AGN on the basis of compact 
jet morphology; background AGN on the basis of 
compact morphology; background AGN on the ba- 
sis of compactness and position. 

3.2.1. Background AGN on the basis of large- 
scale morphology 

37W158 and 37W093 appear to be the cores of 
double-lobed radio sources which are clearly visi- 
ble in all low- resolution maps. We no te, however, 
that B90 and IWalterbos et all (|l985[) treated the 
core and the lobes of 37W093 as separate sources, 
and that the source is located close to the star- 
forming ring. 37W093 also seems likely to be as- 
sociated with the X-ray source Stiele 715, though 
curiously it is offset by 18" from the VLBI position 
(a 2.8(7 error on the X-ray position). 

3.2.2. Background AGN on the basis of compact 
jet morphology 

Three of our sources (37W144, 37W123 and 
37W115) show a clear core jet morphology in the 
VLBI maps and we classify them as background 
AGN on that basis. 



derived position and morphology of the source are 
entirely consistent with our observations, however 
the flux is significantly lower. 

Gelfand et all (|2005l) suggest this source as a 



Pulsar Wind Nebula (PWN) candidate on the ba- 
sis that it is within the optical disk and has a 
supernova remnant (SNR) classification (though 
the candidate SNR optical counterpart is almost 
1' away). The fact that most of the B90 flux is 
concentrated within 4 mas (=0.014 pc) strongly ar- 
gues against this being a Pulsar Wind Nebula. 

also classiiy this source as a 



Gelfand et al 



High-frequency variable, since it v aries by (21 zfc 
2)% between the B90 and NVSS (|Condon et al.l 
Il998l ) surveys. Its 325 MHz flux also varies by 
(51 ± 15)% between the GLG and WENSS sur- 
veys. We see further evidence for this variability 
in that the current VLA 3.4 GHz flux is 30% higher 
than the Braun flux. There is also a difference in 
brightness of a factor of two between our VLBA 
data and that from the VLBA archive. This level 
of variability is consistent with the source being a 
background AGN. 

37W123 

37W123 shows a core-jet structure with the jet 
some 8 mas east and slightly north of the core. 



37W144 

37W144 has one of the highest flux densities of any 
of our sources and was used as an in-beam calibra- 
tor. Though it appears to be only weakly resolved 
in our observations, we we do not recover all of the 
flux seen in the B90 catalog. We also had access to 
VLA images at multiple frequencies and data from 
recent 8.4 GHz VLBA observations of a micro- 



quasar within M31 (iMiddleton et al.l 120131 ). The 
source is entirely unresolved in the VLA images 
and is almost flat spectrum (integrated fluxes are 
33.2, 28.1, and 28.1 mJy at 3.4, 5.3, and 7.5 GHz 
respectively). 

Self-calibrated 8.4 GHz VLBA images revealed 
a bright core with a flux density of 10 mJy beam~^ 
with an extension to the Sou th some 4 mas in ex- 
tent (|Middleton et al.l l2013l Supp. Fig. 1). A 
small amount of additional data on this source was 
available in the VLBA archive (see Table |4]) . The 



37W115 

This source is unresolved in B90, however 2 cm 
VLA maps (0039-f412 in iRossetti et"alll2006i as 
one of the B3-VLA CSS sample) reveal a core 
with evidence of a jet extension and two lobes. 
All of these features (core 5.22±0.52mJy, jet 
3.83±0.38mJy, northern lobe 23±2.3mJy, and 
southern lobe 0.84±0.08mJy) are visible in our 
map (Figure [5]) , however we recover only 12% of 
the B90 flux. It is very likely that most of this 
missing flux in the lobes which are largely resolved 
out by the VLBA, leaving only small hot-spots 
visible. 

Gelfand et al.l classify this source as a "high- 
frequency variable" though it only varies by (6 ± 
4)%. Since the lobes would not be expected to be 
variable, this would represent an extremely high 
level of variability were it attributable only to the 
core. It is more likely to be due to experimental 



14 



3.2.3. Background AGN on the basis of compact 
morphology 

37W142 and 37W125 are both located close to 
the core of M31 (1.1' and 4.8' from M31* respec- 
tively), while 37W150 is located on the edge of one 
of the largest H il regions in M31. In addition, 
the diameter of the circular Gaussian fit of these 
three sources is significantly larger than those of 
the other sources (though the fit for 37W150 is 
extremely poor, as evidenced by the large error). 

37W142 

The visibilities for both our data and the archival 
data fit a circular Gaussian. In the images there 
is some suggestion of extended structure, however 
none exceeds 4cr and it is nonetheless confined 
within some 30 mas. 

The VLBA flux matches the B90 flux, thus if 
this were a hot-spot on the edge of a SNR, there is 
no further radio emission. The lack of any X-ray 
counterpart also argues against this hypothesis. 



37W125 is classified by iGelfand et al.l as both a 
"steep spectrum" source (a = —0.45 ± 0.34 and a 
"normal galaxy" based on its 2MASS counterpart 
and moderate spectral curvature {(p — 0.13±0.05). 
In our data it is well-described by a circular Gaus- 
sian 10 mas in extent. The inferred spectral index 
and curvature are small compared to their errors 
and could be due to source variability. 

37W150 

This source lies just on the edge of one of the 
larger and brighter H ii regions in the ring. It lies 
extremely close to two H ii regions cataloged by 
Walterbos fc BraunI (|1992l ). Although the source 
is only weakly detected, there is evidence of some 
extended structure, possibly jet-like. 

It has an X-ray counterpart which is classified 
as an AGN based on its X-ray spectrum. 

3.2.4. Background AGN on the basis of compact- 
ness and position 

The remaining sources (37W156, B157, 37W092, 
B23, 37W129, 37W080, 37W118, 37W060) are un- 



resolved and are located well away from the core, 
starforming ring and significant H 11 regions of 
M31. Several of them show some evidence of jet 
morphology. These are all likely to be background 
AGN. 

5 of these sources have fluxes (as measured by 
BLOBCAT) significantly (> 1.5a) brighter than that 
in B90. Two of these five sources have previously 
shown evidence of variability. 37W118 has a 37W 
flux which is 2.7(7 higher than the B90 flux but 
consistent with our own. Likewise, 37W129 has 
an NVSS flux which is 1.7cr higher than the B90 
flux but consistent with our own. 

This is a much higher fraction of sources 
with VLBI flux > VLA flux than that found by 
Middelberg et all ()2013i see Figure 14). This may 
be due to the two-decade gap between the B90 ob- 
servations and our own. An alternate possibility 
is that the error on the absolute flux density scale 
in one or the other catalog is underestimated. 

3.3. Stability of ICRF J003824.8+413706 

Figure [S] shows images of the phase reference 
calibrator in archival data compared with our 
data. A modest increase in flux is seen between 
epochs (301±30 mJy vs. 386±39 mJy) and there 
is some evidence for evolution of source structure. 

This is reflected in the change in the ICRF posi- 
tion between the epochs of approximately 0.3mas 
SE. In contrast, we have several sources which 
show no extended structure whose current mea- 
surement errors are significantly smaller than this. 
The use of these sources, or an ensemble of them 
is likely to produce a much more robust reference 
frame for multi-epoch astrometric observations of 
M31 than the ICRF source alone. 

3.4. Possible scatter-broadening by the 
ISM of M31 

Figure [7] shows the size of each source as deter- 
mined by model fitting of the visibilities, plotted 
against its projected distance from M31* (37W123 
and 37W115 are excluded from this analysis on 
account of their mas-scale extended structure. 
37W150 is excluded on the basis of the poor con- 
straint of its fit). It is immediately obvious that 
the two sources closest to M31* (37W142 and 
37W125) have the largest size. 

Images of 37W142 and 37W125 show no ev- 
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(a) (b) (c) 



Fig. 6.— A comparison of the phase reference source ICRF J003824.8+413706 at 18 cm for (a) VLBA 
archival data (Epoch 2003-September-19) with contours starting at ±3x the rms value of 0.47 mJy; (b) our 
data (Epoch 2010-07-4) with contours chosen to match the former; (c) our data with contours starting at 
±3x the rms value of 0.089 mJy. Contours progress in factors of 2. 



idence of complex structure (such as the jet of 
37W123). On the contrary, there is no evidence 
for higher spatial frequencies in the source bright- 
ness distribution at all. 

One possible explanation for this relation is 
that the sources are being scatter-broadened by 
turbulence in the ionized interstellar medium 
(ISM) of M31. Essentially, each of our visibili- 
ties is a direct measurement of the correlation in 
the wave field on transverse separations of between 
100km and 8600km (our UV range). It therefore 
measures the correlation in the plasma column 
density through M31 on these scales — equivalent 
to an angular resolution of 1-100 picoarcseconds. 

To determine whether this is a plausible ex- 
planation, we t ake the main compon ents of the 
NE2001 model (|Cordes fc Lazioll2002l) . which de- 
scribes the ionized ISM of our Galaxy, and use this 
to predict, to an order of magnitude, the scat- 
tering within the inner part of M31. Following 
Tavlor fc CordesI (|1993I ) , the scatter-broadened 
angle (in mas) for a source external to our Galaxy 
is given by 



128 SM°-^ V- 



2.2 



(7) 



where v is the frequency in GHz and SM ^ the scat- 
tering measure is given by 



SM= TleiSfFdS, 



(8) 



where S is the line of sight and F is a scaling 
factor which relates the turbulence to the electron 
density squared. 

3.4-1- Scattering due to a galactic center compo- 
nent 

Observations show the existence of an ex- 
tremely turbulent region close to Sgr A*. In 
NE2001 this region has a height and radius of 
^GC=0.145kpc and iJGc=0.026kpc respectively. 
Both the Ue and F are orders of magnitude higher 
than elsewhere in the galaxy. 

At 1.6 GHz, the scatter-broadened size of a 
background source seen through this region is of 
order 1'. Therefore the detection of 37W142, at 
a projected distance only 0.25 kpc from the core 
already places a constraint on the extent of such 
an extreme scattering region around M31*. 

This result i s cons istent with the findings of 
CiarduUo et al.l (|l988l ). who probed the electron 
density of the inner region of M31 via S ii line ra- 
tios. They found a drop of several orders of mag- 
nitude in electron density over a range to 0.2 kpc 
in radius. 

3-4-2- Scattering due to a 'thin' disk 

Outside the core, the large-scale component of 
NE2001 with the highest scattering measure is 
the inner 'thin' disk (with the 'thick' disk and 
spiral arms being more diffuse and with much 
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Fig. 7. — Size of sources as determined by circular Gaussian fit plotted against distance from M31*. 37W123, 
37W115 and 37W150 are excluded from this plot due to their extended structure (37W150 is also excluded 
due to the large error on its size). 



lower turbulence). This disk is given an annu- 
lar morphology, originally to match a feature seen 
in CO, but favored over a morphology derived 
from the distribution o f ultra-compact H ii regions 
(|Cordes fc Lazioll2003l) . 

The electron density of this component varies 
by two orders of magnitude over the inner 5 kpc of 
the Galaxy. However taking the peak characteris- 
tic density we find that it could give rise to scatter 
broadening of '--^SOmas at 1.6 GHz. The observed 
source sizes that we see are therefore consistent 
with a scattering measure similar to that inferred 
for the inner regions of our own galaxy. Further- 
more, it is possible that some of the scatter that 
we see for sources further than 5' from M31* is 



also due to scatter broadening. 

While the circular Gaussian fit is very poorly 
constrained due to the low S/N, the source 
37W150 does appear to be resolved (Figure [3]). 
If this is due to scatter broadening, rather than 
intrinsic source structure, it would be unsurpris- 
ing given this source's proximity to a large H ii 
region. 

4. Discussion 

4.1. Fraction of sources detected 

As discussed in Section[31 16 individual sources 
were detected (see Figure [5]). If we consider any 
field within 2.5" of another a duplicate, this gives 
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us a detection rate of 16/217. Recall, however, 
that no flux limit was included in our source se- 
lection criteria, for reasons set out in Section [23] 

For more realistic estimate detection we restrict 
ourselves to the B90 catalog and consider only 
those sources that would be detected at 6.6(7 or 
more (assuming no flux is resolved out). This gives 
us a detection rate of 14/97. 

Finally by considering sources that were unre- 
solved or only very weakly resolved (major axis 
<6"), we find a detection rate of 13/48. This 
is broadly consis t ent w ith the detection rate of 
Middelberg et al.l (|2013[) . 

The 35/48 sources in this sample that were not 
detected are also shown in Figure HI There appear 
to be two double sources, a further three are asso- 
ciated with compact H ii regions and are therefore 
likely to be intrinsic to M31 and many appear to 
be clustered in a region close to the starforming 
ring. 

Non-detection of a large number of objects in 
the starforming ring unresolved by the VLA in- 
dicate that a large number of sources intrinsic to 
M31 could be imaged by instruments intermedi- 
ate in resolution between the VLA and the VLBA 
(such as long-baseline LOFAR, e-MERLIN and 
the SKA). 

4.2. Future Work 

The Galactic ionized interstellar medium and 
its turbulence are a subject of intense study. Tech- 
niques used to probe the interstellar medium and 
the ap plications of an accurate model are enumer- 
ated in'Cord es fc Laziol (|2002l ). 

The difficulty of mapping the electron density 
distribution in different parts of the Galaxy from 
our position well within the plane should be noted. 
For example, as recently as 1998 there was signif- 
icant uncertainty over whether the region of en- 
hanced scattering responsible for the observed an- 
gular size of Sgr A* was at the Galactic center 
or somewhere in the plane along the line of sight 



(jLazio fc Cordeslll99S 

Utilizing the recent VLBA sensitivity upgrade 
and the techniques outlined here, it would be pos- 
sible to probe in the ionized medium of another 
galaxy, in detail, for the first time. Furthermore, 
scatter broadening can be discerned from intrin- 
sic size (which scales with A to the power of one 



or less) by its dependence, allowing scatter- 
broadening to be quantified even in cases where 
there is small-scale structure. 

We have scheduled multi-frequency wide-field 
VLBI observations with the aim of finding unam- 
biguous evidence that the observed source sizes 
are due to scatter broadening; increasing the S/N 
on our detections; detecting a larger population 
of sources; and discriminating between scatter 
broadening and source structure for all sources. 
This will also allow us to explore the apparent 
variability of some of the sources. 

As the next generation of radio telescopes 
come online, scatter-broadening of background 
sources will become a remarkably powerful probe 
of the ionized ISM o f other galaxies. Simulations 
(jWilman et al.ll2008l) predict that at 1.4 GHz there 
are approximately 2000 compact sources brighter 
than 10|J.Jy per square degree. Thus with only an 
order of magnitude sensitivity improvement it will 
be possible to probe the ISM of M31 and other 
nearby galaxies in exquisite detail. Furthermore, 
the long-baseline SKA will also operate down to 
500MHz meaning that it will be sensitive to scat- 
tering measures an order of magnitude weaker 
than those probed here. 

The absolute positions of the 16 detected 
sources in M31 have been established with mas 
accuracy. Relative to the in-beam calibrator 
source 37W144, however, the positional accu- 
racy is much higher, in the range of several tens 
to several hundreds of microarcseconds, depend- 
ing on angular separation and source brightness. 
For most sources the S /N-limited fitting accuracy 
dominates the error budget, rather than system- 
atic differential calibration errors. In addition 
to 37W144, the sources 37W115, 37W158 and 
37W092 are all bright and compact enough to 
serve as VLBI calibrators, and these 4 sources 
together could form a calibrator grid capable of 
supplying a relative reference frame accurate at 
the level of 10s of |i.as for future astrometric stud- 
ies of weak sources intrinsic to M31. 

4.3. Conclusions 

• We have detected 16 compact sources across 
a region 35' in diameter. 

• We have classified all 16 as background AGN 
based on their high compactness and com- 
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pact and extended jet morphologies. 

• The relative source positions that we have 
derived (with an accuracy '-^ 10 pas) can be 
used as the basis of a reference frame for any 
future astrometric studies of M31. 

• Our source detection rate is broadly con- 
sistent with other wide-field VLBI studies. 
This rules out a large population of back- 
ground objects scattered, absorbed or oth- 
erwise obscured by the intervening galaxy. 

• We have probed the interstellar medium of 
another galaxy (M31), along several lines of 
sight, for the first time, and have found ev- 
idence of scatter broadening along the two 
lines of sight which pass closest to M31*. 

• Upper limits on scatter broadening along the 
other lines of sight are broadly consistent 
with models for the ionized ISM of our own 
Galaxy. 

• The detection of the source 37W142 places 
an upper limit of 0.25 kpc on the radius of 
any extreme scattering region surrounding 
M31*. 

• Measurements of scatter broadening of back- 
ground sources at multiple GHz frequencies 
with wide-field VLBI has potential as a pow- 
erful probe of the ionized ISM of nearby 
galaxies. 

We would like to thank Walter Brisken (NRAO) 
for assistance in data transfer. We thank Stefanie 
Muehle and Bob Campbell (JIVE) for useful dis- 
cussions regarding VLBI scheduling and correla- 
tion. Extensive use was made of topcat H ( Tavloi 
20051 ) and Aladiifl (|Bonnarel et al.ll2000[ ). The 
National Radio Astronomy Observatory is a facil- 
ity of the National Science Foundation operated 
under cooperative agreement by Associated Uni- 
versities, Inc. 
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Fig. 5. — Image of 37W115 with a uv cutoff of 
lOMA to show large-scale structure. 



This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.2. 
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Fig. 8. — Grayscale: image extracted from the NVSS survey (|Condon et al.lll998l ). Red 'x' with blue labels: 
detected sources. Red '+': unresolved B90 sources (major a xis > 6") that would be de tected at 6.6 u were 
all their flux com pact. Green di amonds: H ii regions from IWalterbos fc BraunI (jl992l ). Pink circles: H ii 
regions from .Azimlu et al.l (|2011l ). Black square: pointing center. 



22 



